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LONG  TERM  GOALS 

The  long  term  goal  of  this  project  is  to  develop  efficient  inversion  algorithms  for  successful 
geoacoustic  parameter  estimation,  inversion  for  sound-speed  in  the  water-column,  and  source 
localization,  exploiting  (fully  or  partially)  the  physics  of  the  propagation  medium.  Algorithms  are 
designed  for  inversion  via  the  extraction  features  of  the  acoustic  field  and  optimization.  The  potential 
of  analytic  approaches  are  also  investigated. 

OBJECTIVES 

•  Achieve  accurate  and  computationally  efficient  inversion  for  propagation  medium  parameters  and 
source  localization  by  designing  estimation  schemes  that  combine  acoustic  field  and  statistical 
modeling. 

•  Develop  methods  for  passive  localization  and  inversion  of  environmental  parameters  that  select 
features  of  propagation  that  are  essential  to  model  for  accurate  inversion. 

•  Implement  Bayesian  filtering  methods  that  provide  dynamic  and  efficient  solutions  for  the  first  two 
objectives. 

•  Develop  analytic  techniques  for  sediment  sound  speed  estimation. 

APPROACH 

Continuing  efforts  from  previous  years,  we  worked  with  Bayesian  approaches  applied  to  sound  signals 
for  the  extraction  of  acoustic  features  using  a  combination  of  physics  and  statistical  signal  processing. 

One  of  the  topics  approached  this  past  year  was  source  localization,  bathymetry,  and  water  column 
sound  speed  estimation  using  arrival  time  estimates  for  propagation  in  multipath  environments  with 
sequential  Monte  Carlo  methods,  namely,  particle  filtering,  tied  with  a  linearization  method  with  novel 
features.  Some  of  our  results  have  been  presented  in  [1,  2,  3,  4]  and  [5];  the  latter  paper,  focusing  on 
the  linearization  approach,  is  currently  under  review.  The  method  is  an  extension  of  previous  methods 
developed  by  the  PI,  now  using  more  advanced  techniques.  The  initial  goal  is  to  estimate  accurately 
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arrival  times  of  sound  paths  in  shallow  water  environments.  Then,  we  propagate  these  arrival  times 
and  their  posterior  probability  density  functions  (PDFs)  through  a  quasilinear  model  for  source 
location,  bathymetry,  and  water  column  sound  speed  profile  estimation. 

More  work  with  a  new  statistical  model  more  accurately  describing  dispersion  in  the  ocean  was 
conducted  for  identifying  arrival  times  and  amplitudes  of  distinct  frequencies  (within  a  single  mode  or 
across  different  modes).  Those  characteristics  provide  significant  information  on  properties  of  the 
propagation  medium  and  source  location.  The  importance  of  modal  arrival  times  and  amplitudes  in 
geoacoustic  inversion  and  source  localization  using  dispersion  curves  has  been  extensively  discussed 
in  [6,  7,  8,  9]. 

Finally,  we  worked  on  a  new  sediment  sound  speed  estimation  scheme  based  on  Stickler’s  inverse 
problem  approach  [10,  1 1].  As  we  will  mention  in  the  next  section,  we  managed  to  estimate  sediment 
sound  speed  in  a  synthetic  waveguide  using  an  analytic  method,  improving  on  previously  obtained 
results  with  a  similar  but  simpler  approach. 

WORK  COMPLETED 

Previous,  approaches  developed  by  the  PI  were  improved  to  provide  more  accurate  estimates  of  arrival 
times  and,  subsequently,  of  source  location,  bathymetry,  and  the  water-column  sound  speed  profile. 
Specifically,  in  addition  to  conventional,  forward  particle  filters,  smoothing  backward  mechanisms 
were  employed,  that  reduced  uncertainty.  The  arrival  time  tracking  method  was  applied  to  SW06  l fin 
data.  It  was  combined  with  a  new  linearization  technique  that  now  includes  differentiation  with 
respect  to  EOF  coefficients  rather  than  sound  speed  directly.  Source  localization,  bathymetry,  and 
sound  speed  profile  results  were  obtained  in  terms  of  both  point  estimates  and  posterior  PDFs,  which 
will  be  shown  later. 

Additionally,  the  Pi’s  particle  filtering  method  was  applied  to  both  synthetic  and  real  long-range  data 
for  dispersion  curve  extraction  with  a  new  and  more  accurate  statistical  model  than  the  simple 
Gaussian  one  that  had  been  previously  used.  Bathymetry  estimates  were  extracted  using  the  estimated 
dispersion  curves  for  the  synthetic  case  (more  work  is  currently  conducted  with  the  real  data). 

Finally,  we  estimated  sediment  sound  speed  using  a  new  analytic  method,  building  on  work  presented 
in  [10,  11]. 

RESULTS 

As  mentioned  in  the  previous  section,  we  used  the  particle  filtering  method  along  with  a  linearization 
scheme  that  related  arrival  time  estimates  to  source  range  and  depth,  water  column  depth,  sound  speed 
profile  in  the  water  column,  and  tilt.  The  Jacobian  matrix  necessary  for  linearization  can  become 
complex  if  differentiation  with  respect  to  the  depth  dependent  sound  speed  is  perfonned.  We 
developed  an  alternative  scheme,  where  the  Jacobian  matrix  includes  derivatives  with  respect  to  EOFs 
instead  of  sound  speed  directly.  PDFs  for  source  localization,  water  column  depth,  and  sound  speed 
profile  EOF  coefficients  obtained  with  this  method  are  shown  in  Figure  1 .  The  method  required  only  a 
few  iterations  for  convergence  and  the  results  were  in  agreement  with  estimates  obtained  at  nearby 
sites  [12,  13]  with  more  computationally  demanding  global  optimization  techniques.  It  should  be 
noted  that  no  inversion  results  have  been  previously  presented  from  the  Ifm  data  used  in  this  work. 
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Figure  2  shows  (a)  the  mean  sound  speed  profile  as  calculated  by  CTD  measurements  at  the  site  and 
the  sound  speed  profile  constructed  with  the  MAP  estimates  obtained  from  the  extracted  arrival  time 
PDFs  and  linearization  and  (b)  the  profiles  of  (a)  with  superimposed  sound  speed  profiles  constructed 
from  distinct  arrival  time  particles. 

Figure  3  shows  the  real  data  time-series  (solid  lines)  and  synthetic  time-series  (dotted  lines)  generated 
using  the  parameter  estimates  from  the  filtering/linearization  scheme.  The  synthetic  time-series  agree 
closely  with  the  true  data,  strongly  indicating  that  our  estimates  are  close  to  the  true  values. 

It  was  previously  mentioned  that  work  was  perfonned  on  extracting  dispersion  curves  for  long-range 
propagation  with  a  new  statistical  model.  Preliminary  results  show  that  the  new  model  offers  higher 
reliability  in  the  estimation  of  dispersion. 

Lastly,  Figure  4  shows  sound  speed  estimation  in  sediments  employing  the  method  mentioned  in  the 
past  section.  The  interpolant  approach  developed  in  our  work  provides  closer  profiles  to  the  true  one 
than  other  approaches,  with  the  improvement  mostly  evident  at  the  interfaces  between  different 
sediments. 
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Figure  1:  PDFs  of  (a)  source  range,  (b)  source  depth,  (c)  water  column  depth, 

(d)  pi,  (e)  p2,  (f)  /<  ;  for  real  data. 


3 


0 


a.  40 

03 

~G 

5  60 


1480 


1490 


1500 


1510  1520 

Sound  speed  (m/s) 


1530 


1540 


(a) 


1550 


Sound  speed  (m/s) 


Figure  2:  (a)  The  mean  sound  speed  profile  as  calculated  by  the  CTD  measurements  (squares)  and 
the  sound  speed  profile  constructed  with  the  MAP  estimates  obtained  from  the  extracted  arrival  time 
PDFs  (circles),  (b)  The  profiles  of  (a)  with  superimposed  sound  speed  profiles  constructed  from 

distinct  arrival  time  particles  (solid  curves). 
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Figure  3:  Real  data  (solid  lines)  and  synthetic  time-series  (dotted  lines)  generated  using  the 

linearization  estimates. 
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Figure  4:  Sediment  sound  speed  estimation  with  an  analytic  method,  improving  on  results 

previously  obtained  from  similar  techniques. 


IMPACT 

Accurate  arrival  time  estimation  and  fast  inversion  methods  for  source  localization  and  water-column 
sound  speed  estimation  have  been  studied  with  the  goal  of  producing  accurate  parameter  estimates  and 
measures  of  the  uncertainty  in  the  inversion  process,  new  approaches  facilitate  arrival  time  extraction 
and  the  association  between  paths/modes  and  detected  arrivals  and  also  produce  PDFs  of  arrival  times, 
modal  frequencies,  and  corresponding  amplitudes.  These  are  critical  features  in  estimating  inversion 
parameters  efficiently  and  accurately.  The  linearization  method  also  impacts  favorably  inversion 
adding  to  the  efficiency  and  accuracy  just  mentioned.  Also,  the  development  of  a  new  inversion 
technique  has  a  significant  impact  on  the  problem  of  geoacoustic  parameter  estimation,  because  of  its 
potential  to  estimate  sound  speed  in  sediments  in  an  exact,  fast,  analytic  fashion. 

RELATED  PROJECTS 

The  PI  is  collaborating  with  Drs.  Yardim  and  Gerstoft  on  sequential  filtering  in  ocean  acoustics.  The 
PI  is  also  collaborating  with  Dr.  Leon  Cohen  on  comparing  numerical  and  analytical  descriptions  of 
dispersion  in  ocean  environments  and  is  involved  in  discussions  with  Dr.  Ross  Chapman  on  inversion 
for  attenuation. 
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